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Detection  of  nerrow- bend  low-level  signals! 

Simulated  performance  of 
high-resolution  eigenstructure  methods 
versus  conventional  beamforming 

S.  Jesus 


Executive  Summary!  In  many  instances,  sonar  performance  and  hence 
ASW  capability  is  limited  by  the  sonar’s  spatial  resolution,  i.e.  its  ability 
to  separate  aconstic  signals  arriving  from  close  angles  in  order  to  obtain 
both  noise  rejection  and  angle  of  arrival  for  the  location  and  classification 
of  targets.  This  limitation  is  directly  related  to  the  length  of  the  array,  the 
aperture. 

In  the  last  two  decades  a  number  of  signal  processing  techniques  have  been 
proposed  that  increase  the  resolution  power  of  a  given  array  without  increas¬ 
ing  its  length:  these  are  the  so-called  high-resolntion  techniques.  However, 
although  successful  in  other  fields,  their  adoption  in  operational  sonar  sys¬ 
tems  has  been  slow.  One  reason  for  this  is  the  difficulty  in  determining 
the  detection  capability  of  high-resolntion  techniques  compared  with  con¬ 
ventional  techniques,  due  to  the  absence  of  a  commonly  accepted  measure 
of  performance  that  can  be  expressed  in  operational  terms. 

The  aim  of  the  present  work  is  to  improve  this  situation  by  showing,  with 
a  simulation  study,  that  the  use  of  an  appropriate  high-resolution  technique 
may  give  a  higher,  and  sometimes  much  higher,  detection  performance  than 
the  conventional  beamformer.  Improvement  can  be  obtained  especially  when 
detecting  a  low-level  target  in  the  vicinity  of  one  or  more  high-level  interfer¬ 
ing  targets.  This  situation  can  occur  in  many  crucial  geographic  areas  where 
the  acoustic  noise  field  is  in  major  part  due  to  the  contribution  of  individual 
ships  behaving  as  interfering  targets  generating  ship-indnced  noise  fields. 

The  results  of  this  work  contribute  to  the  understanding  of  the  detection 
behaviour  of  the  high- resolution  techniques  in  a  number  of  tealistic  field 
situations.  Such  techniques  can  also  be  viewed  as  an  alternative  to  increas¬ 
ing  the  array  length  and  thereby  improving  the  operational  characteristics 
of  existing  sonar  systems.  A  companion  report  in  preparation  shows  the 
detection  performance  of  these  high-resolution  techniques  using  real  data. 

Further  progress  requires  improved  knowledge  of  the  background  noise  field, 
its  parametrisatlon  in  a  noise  model  and  its  inclusion  is  the  noise  source  es¬ 
timation  process  as  a  priori  information.  The  high-resolution  techniques 
presented  here  allow  such  information  to  be  included.  The  results  may  sub¬ 
stantially  improve  if  this  information  is  correct. 
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V 

Abstract:  In  this  report  the  detection  of  narrow-band  low-level  signals 
by  passive  sonar  is  addressed,  specifically  the  detection  of  a  low-level  point 
source  in  a  noise  field  modelled  as  a  large  number  of  high-level  source  lines 
embedded  in  white  additive  noise.  This  model  may  give,  with  a  conven¬ 
tional  processor,  a  measured  noise  field  where  the  low-level  source  can  be 
masked  by  high-level  sources  along  similar  directions.  In  this  context,  high- 
resolution  eigenstructure  detection  algorithms  are  shown  to  achieve  a  signif¬ 
icantly  higher  probability  of  detection  than  the  conventional  beamformer. 
The  performance  of  four  detection  algorithms  is  compared  using  synthetic 
data  -  a  companion  report  will  present  real  data  results. 
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DOA  (dag) 


Fig.  1.  Estimated  power  (dB)  vs  bearing  (deg)  -  conventional 
beamformer  with  Hann  shading  (a)  and  high-resolution  eigenstruc- 
ture  method  (b):  three  sources  at  locations  -10°,  -5°  and  0°; 
respective  powers  30  dB,  0  dB  and  30  dB;  number  of  averages 
N  =  200,  16  sensors. 


Simulated  tests  have  been  conducted  with  different  noise  fields  for  variable  relative 
location  of  the  low-level  source  to  the  high-level  interfering  source(s),  for  variable 
signal-to-noise  ratios,  for  different  array  shadings  (conventional  beamformer)  and 
also  for  different  observation  times  (‘time-bandwidth  products’).  In  a  companion 
report  [7]  the  methods  described  here  are  applied  to  real  data  and  compared  with 
previous  results  [6]. 
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2.  Background 

Passive  sonar  target  detection  and  classification  is  generally  carried  out  by  identify¬ 
ing  the  spectral  content  of  a  given  data  set  with  known  ship  and  submarine  acoustic 
spectral  signatures.  Source  location  is  achieved  by  spatial  processing  to  extract  the 
directional  features  of  the  data.  Time  and  space  operations  are  generally  performed 
separately.  In  this  study  we  will  concentrate  on  the  spatial  data  processing  for  source 
detection  and  direction-of-arrival  (DOA)  estimation.  Signals  of  interest  are  assumed 
to  be  of  the  narrow-band  type  for  which  the  frequency  bandwidth  A/  is  less  than 
or  equal  to  the  frequency  resolution  1/To  of  the  fast  fourier  transform  (FFT)  used 
for  spectrum  estimation,  A /  <  1/To. 


2.1.  DATA  MODEL 

The  spatial  field  is  sampled  by  a  line  array  of  L  omnidirectional  sensors  at  a  constant 
interval  d.  To  avoid  spatial  aliasing,  d  is  assumed  to  be  always  <  c/2/m„,  where  c 
is  the  sound  speed  and  /,„»*  is  the  highest  frequency  contained  in  the  signal.  Let 
*/(t'T,)  be  the  /th  sensor  output  at  a  discrete  time  iT,  where  the  sampling  period  T, 
is  chosen  such  that  T,  <  l/2/m„  The  signal  i;(iT,)  is  assumed  to  be  stationary 
over  a  record  length  T,.  This  interval  T,  is  then  segmented  into  N  anal  intervals 
with  2 K  samples  each,  T,  =  2N KT,.  Each  of  the  intervals  is  then  passed  through 
an  FFT.  Let 

yT(n,  /ic)  =  [»>(",/),),  y2(n,/fc),.  .,yL(n,fk)},  n  =  1 - AT;  fc  =  1, . . . ,  K , 

(1) 

be  the  L-dimensional  array  complex  vector  at  time  ‘snapshot’  n  and  for  frequency 
fk\  T  stands  for  the  transpose.  N  is  termed  the  ‘time-bandwidth  product’,  which 
is  a  compromise  between  data  stationarity  and  power  spectrum  estimate  stability. 
The  spatial  processing  of  y(n,  fk)  will  be  achieved  separately  for  each  frequency  fk 
Keeping  this  in  mind,  the  frequency  index  will  be  omitted  in  the  following  for  clarity. 

The  standard  model  to  represent  the  observation  vector  y(n)  with  I  far-field  sources 

is  [4] 

y(n)  =  As(n)  +  v(«).  (2) 

Matrix  A  is  an  L  x  /  complex  matrix,  the  columns  of  which  are  the  mode  vectors 

aT(0,)  =  [l,e '*’'•* _ ,  , _ /  (3) 

where  the  ith  source  wavenumber  Vqx  is  defined  bv 

2  ir 

Vg_  -  —  sin  6, 
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with  A  being  the  source  wavelength  and  0t  the  angle  defined  between  the  ith  source 
wavefront  and  the  normal  to  the  array.  In  this  case 

A  =  [a(fli  ),a(0j a(0/)].  (4) 


In  the  standard  model  (2),  the  vector 


s(n)  =  (si(n),  j2(n), . .  .  ,S/(n)J 


(5) 


is  the  /-dimensional  signal  complex  (analytic)  vector  wheie 

«i(n)  =  a,t(n)  +  >sfl)(n),  i  =  l,.  ,/.  (5a) 

The  real  and  imaginary  parts  of  si(n)  are  assumed  to  be  uncorrcl^ted,  stationary, 
and  ergodic  gaussian  random  processes  with  zero-mean  and  variance  .  The 
covariance  matrix  of  the  signals  is  denoted 

S  _  £(s(-)s(-)H],  (5b) 

where  H  stands  for  the  conjugate  transpose. 

Finally,  the  complex  additive  noise  v(n)  is  assumed  to  be  stationary,  ergodic,  zero- 
mean  and  with  covariance 

£[v(  )v(  )H|  =  <t,2I.  (6) 

where  I  is  the  identity  matrix  and  is  a  scalar  representing  the  noise  power. 
Moreover,  the  signals  are  assumed  to  be  uncorrelated  with  the  noise  process. 


In  a  DOA-estimation-only  system,  one  uses  the  spectrum  analyzer  output  data 
(y(n);  n  =  1, .  . .  /V }  to  estimate  the  set  of  angles  {0;;i  =  1,...,/}.  The  number 
of  sources  /  is  in  general  not  known  and  also  needs  to  be  estimated.  This  problem, 
when  the  number  of  sources  is  not  known,  is  often  called  a  detection  problem. 
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2.2.  CONVENTIONAL  BEAMFORMER 

Only  the  case  of  a  linear  equispaced  array  will  be  treated.  Conventional  beamforming 
is  a  common  nonparametric  and  computational  efficient  technique  to  combine  the 
sensor  outputs  in  this  particular  array  arrangement.  The  computational  efficiency  of 
the  conventional  beamformer  is  mainly  due  to  the  use  of  a  spatial,  fourier  transform 
implemented  by  the  FFT  algorithm  [8] . 

The  beamformer  operates  on  the  array  by  coherent  summation  of  the  L  sensor 
outputs  with  the  apropriate  spatial  delay  to  steer  a  ‘beam’  to  a  given  direction  9. 
Let 

L 

6n(fl)  =  yi  htwi{e)yi(n)  (7) 

i=i 

be  the  complex  beamformer  output  steered  to  direction  9  at  time  ‘snapshot’  n.  The 
term  contains  the  phase  shift  imposed  to  sensor  i  to  steer  the  beamformer  to 

direction  6, 

w,(9)  =  e  jV  l  1,  (8) 

where  all  quantities  have  been  defined  above.  The  factor  ht  is  a  spatial  weight  asso¬ 
ciated  with  sensor  1;  the  sequence  {hi,  l  =  1  is  known  as  the  spatial  shading 

function.  The  periodogram  power  wavenumber  spectrum  estimate  is  obtained  by 
summing  (7)  over  the  N  time  segments,  which  gives  for  all  beams 

BN=W'HR."W',  (9) 

where  the  matrix  W'  contains  the  direction  vectors  (8)  weighted  by  the  spatial 
function  hi.  The  matrix  R,^  is  a  periodogram  estimate  of  the  data  cross-covariance 
matrix  R.(1) : 

K  -  <10> 

n .-  1 


2.3.  HIGH-RESOLUTION  EIGEN  STRUCTURE  METHOD 

This  method  differs  from  the  previous  one  mainly  in  the  fact  that  it  explicitly  uses 
a  model  of  the  data-generating  mechanism  whereas  the  other  did  not.  It  provides 
much  more  accurate  results  if  the  model  is  appropriate.  This  parametrical  method 
uses  an  interpretation  of  the  measured  data  through  time  as  a  sequence  of  an  L- 
dimensional  vector  moving  in  an  L-dimensional  vector  space.  If  the  array  receives 
a  certain  number  (<  L)  of  coherent  wavefronts  the  movement  of  the  observed  data 
vector  is  restricted  to  a  certain  linear  subspace  of  the  total  L- dimensional  reachable 
space.  The  method  relies  on  the  possibility  of  decomposing  the  total  space  into  two 
linear  subspaces:  the  signal  subspace  and  the  noise  subspace. 
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3.  Source  detection 

Source  detection  is  a  difficult  problem  in  underwater  passive  listening  due  to  the 
interrelation  between  the  errors  introduced  by  the  technique  used  for  array  process¬ 
ing  and  the  test  of  signal  presence  itself.  Assuming  the  data  model  described  in 
Subsect.  2.1  the  optimal  detector  is 

U  pl(y)  =  <2yHPy,  (20) 

where  P  is  the  i-dimensional  complex  matrix 

P  =  (<r,2I  +  ASAh)~‘ASAh.  (21) 


In  practice,  the  detection  characteristics  of  a  receiving  array  are  commonly  described, 
in  statistical  terms,  by  the  probabilities  of  detection  (Pp)  and  of  false  alarm  (Pfa) 
for  a  given  signal-to-noise  ratio  (SNR).  These  two  probabilities  are  respectively  as¬ 
sociated  with  the  detection  of  a  signal  under  the  two  hypotheses 

Hi:  signal  plus  noise  is  present,  or 
Ho:  only  noise  is  present. 

Thus,  the  detection  and  false  alarm  probabilities  are  given  by 

Pp(SNR)  =  Prob{I(y)  >  7(SNR)/H,  },  (22a) 

PFA(SNR)  ^  Prob{/(y)  >  7(SNR)/H0}.  (22b) 

If  !(y)  takes  the  form  (20)  the  detection  processor  is  said  to  be  optimal  (under  the 
assumptions  of  the  data  model  described  in  Subsect.  2  l). 
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3.1.  WITH  THE  CONVENTIONAL  BEAMFORMER 

The  standard  scenario  for  testing  hypothesis  Ho/H,  is  to  consider  the  possibility  of 
presence  of  a  single  target  in  white  additive  noise.  In  this  case,  the  received  power 
estimate  given  hy  the  beamforiner  (9)  directed  at  the  source  location,  say  0, ,  is 

S-ViJlHo  =  *1,  (23a) 

S'ViIIh,  -  +  &l  (23b) 

It  can  be  proved  that  the  conventional  beamformer  output  (23)  reaches,  in  this 
particular  case,  the  optimal  detector  (20).  Both  (23a)  and  (23b)  are  random 
distributed  variables  with  IN  degrees  of  freedom  and  mean  values  crj  and  <r\  + 
a •  respectively.  From  the  respective  densities  it  is  relatively  easy  to  calculate  the 
probability  of  detection  P£  for  a  fixed  threshold  B%.  depending  on  a  given  allowed 
fase  alarm  rate  P ^  (0  stands  for  zero  interfering  sources  -  see  e  g.  (9)). 

Let  us  consider  now  the  more  realistic  scenario  of  the  detection  of  one  low-level  point 
source  in  the  presence  of  a  number  of  high-level  interfering  sources,  all  embedded 

in  white  additive  noise.  Assume  /  point  sources  at  locations  {0,;  i  =  1 . /}.  The 

beamformer  is  steered  to  the  direction  of  the  source  we  want  to  detect,  say  source  1 

at  bearing  0, .  The  /  -  1  interfering  sources  with  powers  {< rj  ;i  =  2 . /}  are 

located  respectively  at  bearings  {0, ;i  =  2,  Clearly  using  (9) 

i 

-  *1  +  £>*.a>"  (24a) 

« =  2 

/ 

#'V,)|h,  -  <  l  *1  +  £<«...  (24b) 

i~2 

where  the  last  term  represents,  in  both  expressions,  the  sum  of  the  'influences’  of  the 
I  1  interfering  sources  on  estimating(and  detecting)  source  1.  These  ‘influences’ 
obviously  depend  on  the  relative  location  of  source  1  to  the  interfering  source(s). 
The  factor  an  is 

1  |  sin  by, (L  |* 

"  L 2 1  sin  bylt  I 

with 

2irfd 

Ti,  -  - (sin0j  -  sm0,). 

r 


The  detector  given  by  (24 )  is  no  longer  optimal  In  effect  ,  with  some  straightforward 
manipulations  it  can  be  shown  that  if  at  least  one  interference  source  is  present, 
i.c.  0l>  /  then  the  probability  of  detection  ip  the  presence  of  interfering 

sources,  Pp,  is  always  less  than  or  equal  to  the  probability  of  detection  with  no 
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interfering  sources,  P&  (for  the  same  false  alarm  probability  PyA  —  /*/, ),  This 
is,  of  course,  due  to  the  beampattern  interference  when  detecting  one  source  in 
a  multi-source  noise  field.  Let  us  recall  that  the  sidelobe  impact  can  be  greatly 
reduced  by  using  well-known  techniques  of  sensor  weighting  (array  shading).  These 
techniques  may  improve  the  probability  of  det.  cting  ‘well-separated’  sources  with  the 
counterpart  of  a  loss  in  resolution  and  therefore  a  loss  of  detecting  ‘closely-separated’ 
sources  (within  the  beampattern  mainlobe). 

In  conclusion,  the  degradation  of  the  detection  characteristics  of  the  conventional 
beamformer  is  due  to  the  array  beampattern  effect  and  is  therefore  related  to  the 
angular  separation  between  sources  and  depends  on  their  relative  amplitude.  We 
will  show  in  Subsect.  3.2  that  the  detection  characteristics  of  the  high-resolution 
eigenstructure  techniques  do  not  suffer  from  any  such  inconveniences:  they  present  a 
constant  detection  behaviour  for  all  directions  and  in  particular  for  low-le'’el  sources 
near  to  strong  interfering  sources. 


3.2.  WITH  THE  HIGH-RESOLUTION  EIGENSTRUCTURE  METHOD 

Making  use  of  the  high-resolution  eigenstructure  method  outlined  in  Subsect.  2.3,  the 
problem  of  source  detection  is  equivalent  to  the  problem  of  estimating  the  dimension 
of  the  signal  subspace  Ei  (17).  The  signal  subspace  dimension  is  equal  to  the  number 
of  highest  non-equal  eigenvalues  of  the  cross- covariance  matrix  Rys/  (12). 

For  a  finite  observation  time,  i.e.  finite  N ,  the  eigenstructure  of  the  estimated  cross¬ 
covariance  matrix  R^,  (10)  based  on  the  observation  set  (1)  is  defined  by 

=  A,e,,  f  -=  1 . L.  (25) 

The  estimated  eigenvalues 

A,  >  A2  >  ...  >  A,  >  ...  >  Al  >  0.  (26) 

are  all  different  with  probability  one  [10]  making  it  difficult  to  find  the  number  of 
signals  by  merely  ‘looking’  at  the  eigenvalue  distribution. 

Two  approaches  have  been  proposed  to  solve  this  problem  -  the  generalized  likeli¬ 
hood  ratio  test  (11)  and  the  information-based  theoretic  criteria  (12): 

•  The  first  approach  splits  the  decision  problem  into  a  sequence  of  yes/no 
answers  to  a  series  of  questions:  ‘Is  the  number  of  sources  equal  to  k  for 
k  =  0,  ...,L  -  1?’  Under  the  hypothesis-testing  theory,  answering  ‘yes’  to 
the  question  is  equivalent  to  accepting  the  hypothesis 

Hj,:  there  are  at  most  k  sources 
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Answering  ‘no’  is  equivalent  to  accepting  Hy: 

Ht  there  are  at  least  k  +  1  sources. 

Starting  with  k  =  0  the  test  is  performed  until  Hy  is  rejected  (equivalent  to 
accepting  H*)  or  fc  =  L  -  1,  in  which  case  k  =  L. 

•  The  second  approach  treats  the  decision  problem  as  the  determination  of  the 
rank  of  the  signal  matrix  which  is  considered  as  a  model  selection  problem. 
The  selected  model  is  the  one  that  ‘best.’  fits  the  data  under  certain  criteria. 
Two  methods  for  this  approach  have  been  proposed:  the  Akaike  information 
criteria  (AIC)  and  the  minimum  descriptor  length  (MDL). 


All  three  methods  referred  to  above  are  ultimately  seeking  to  test  for  the  equality 
of  the  L  -  k  smallest  eigenvalues  of  the  sample  cross-covariance  matrix.  This  test  is 
known  in  the  literature  as  ‘the  sphericity  test’  [13,14]. 


The  assumption  that  the  data  is  a  series  of  zero-mean  statistically  independent 
gaussian  random  vectors  leads  to  an  optimum  estimation  of  the  model  parameters 
provided  by  the  maximum  likelihood  estimator.  It  is  well  known  that  in  this  case 
the  functional 

T(k)=  -21og«t.  (27) 

where  a *  is  the  likelihood  ratio  given  by 


/  Ili=fc+I  Ai _  \ 


(28) 


is  distributed  for  large  samples  approximately  as  \2  (15,16).  This  statistic  is  used  by 
the  three  criteria,  the  difference  being  only  in  the  modifications  used  to  approximate 
the  asymptotical  y2  distribution  for  short  or  moderate  sample  size.  We  are  not 
attempting  to  derive  the  whole  underlying  theory  as  this  has  already  been  extensively 
covered  in  the  literature  referred  to  above;  only  the  expressions  for  the  criteria  are 
given  below. 


Generalized  likelihood  ratio  test  (GLRT)  Following  Liggett 
for  moderate  sample  size  would  be 

(11)  a  corrected  test 

TglrtW  =  l‘(k)  logo,. 

(29) 

with  the  multiplying  factor  /((*■)  equal  to 

...  -2(N  -  k  (2 (L  kf  i  1]) 

"(*)=  6 a  *r 

(30) 

This  test  is  conducted  until 

7'ulrt(*')  >0 

(31) 
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or  k  -s  L  -  1.  The  detection  threshold  is  drawn  from  a  .v*  distribution  with 
i/*  -  (L  -  ib)2  -  1  degrees  of  freedom  for  a  given  significance  rjk  fixed  by  the  operator 
such  that 


ProbU*4  >  *£}  =  Vk- 


(32) 


Akaike  Information  Criteria  (AIC)  From  (12),  the  estimated  number  of  signals  is, 
according  to  the  AIC  criteria,  the  value  of  k  =  0,  . . ,  L  -  1  for  which 

2aic(*)=  -2Aloga),  +  2fc(2L-  it)  (33) 


is  minimum. 

Minimum  Descriptor  Length  (MDL)  Also  from  (12),  the  number  of  signals  is  equal 
to  the  value  of  k  -  0, .  . . ,  L  -  1  that  minimizes  the  functional 

Tmdl(*)  =  -Aloga*+  \k(2L-  k)\ogN.  (34) 


Comparing  these  detection  procedures,  the  first  remark  is  that  they  are  obviously 
independent  of  the  relative  location  of  the  sources,  in  contrast  to  the  conventional 
beamformer.  They  depend  only  on  the  observation  time  or  equivalently  on  the  num¬ 
ber  of  samples  N .  This  is  the  well-known  behaviour  of  the  high-resolution  eigen- 
structure  techniques  which  have  asymptotic  (large  N)  infinite  resolving  power  (5). 
This  will  be  clearly  illustrated  in  Sect  4  by  simulation. 
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4.  Simulation  results 

For  each  case  the  probabilities  of  detection  are  estimated  as  a  function  of  the  input 
signal- to-noise  ratio  (SNRj„)  and  the  relative  location  (DOAr)  of  the  low-level  source 
to  the  high-level  interfering  source)  s).  Two  simulated  noise  fields  were  generated 
to  illustrate  the  comparative  behaviour  of  the  four  detection  algorithms.  In  the 
first  example  the  noise  field  is  formed  by  two  line  sources  in  a  flat  background 
noise:  one  is  the  high-level  interfering  source  which  is  fixed  in  location  and  in  power 
and  the  other  is  the  low-level  source  on  which  the  detection  is  performed.  In  the 
second  example  the  low-level  source  is  located  between  two  high-level  sources,  having 
a  fixed  power  but  a  variable  angular  separation.  The  data  has  been  generated 
according  to  model  (2)  with  the  parameters  given  in  Table  1.  In  both  examples  two 
‘tiine-bandwidth  products’  were  tested.  Also  two  array  shadings  were  used  in  the 
conventional  beamformer:  uniform  and  Hann  shading.  For  each  case  four  curves  are 
plotted  corresponding  to  the  four  detection  algorithms  being  tested: 

•  Conventional  beamformer  (CVB  -  curve  C). 

•  Akaike  information  criteria  (AIC  -  curve  A). 

•  Minimum  descriptor  length  (MDL  -  curve  M). 

«  Generalised  likelihood  ratio  test  (GLRT  -  curve  G). 


4.1.  EXAMPLE  1:  ONE  HIGH-LEVEL  INTERFERING  SOURCE 

The  empirical  detection  characteristics  are  shown  in  Figs.  2  to  5.  The  probability 
of  detection,  for  a  constant  false  alarm  probability,  is  given  vs  the  input  signal-to- 
noise  ratio  in  Fig.  2  and  vs  the  relative  location  in  Figs.  3  to  5.  The  probability  of 
false  alarm  can  be  controlled  in  the  conventional  beamformer  but  not  in  the  other 
detection  algorithms.  In  this  example,  the  empirical  estimation  of  the  false  alarm 
rate  over  000  statistically  independent  draws  gave  the  results  shown  in  Table  2. 
Among  th-  three  high-resolution  detection  algorithms  one  can  note  that  while  the 
AIC  algorithm  provides  the  largest  probability  of  false  alarm  (2.7%),  the  probability 
of  false  alarm  of  the  MDL  algorithm  could  not  be  evaluated  with  this  sample  size. 
The  mean  Pta  of  all  three  high-resolution  algorithms  is  about  1%,  which  is  aprox- 
imately  the  value  for  the  conventional  beamformer;  therefore  the  true  comparable 
characteristics  are  relatively  close  to  the  curves  shown. 


Comments  Observing  Fig.  2,  one  can  remark  that  the  CVB  performs  better  than 
all  the  other  algorithms  for  DOAr  —  -30°  (a),  its  performance  then  gradually 
decreases  when  approaching  the  high-level  source  (b)  and  (c),  and  is  null  for  DOAr  = 
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Fig.  2a.  Estimated  probability  of  detection  (%)  vs  the  input 
signal-to-noise  ratio  (dB),  one  interfering  source  at  -10°,  N  —  50, 
No  —  200,  Hann  weighting  DOAr  -  -30°. 


SNR*  (dB) 


Fig.  2b.  Estimated  probability  of  detection  (%)  t>s  the  input 
signal-to-noise  ratio  (dB),  one  interfering  source  at  10c,  N  =  50, 
No  ~  200,  Hann  weighting  DOAr  =  -  18° 
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Fig.  2c.  Estimated  probability  of  detection  (%)  vs  the  input  signal- 
to-noise  ratio  (dB),  one  interfering  source  at  10°.  N  =  50,  No  = 
200,  Hann  weighting  DOAr  =  -10°. 


Fig.  2d.  Estimated  probability  of  detection  (%)  vs  the  input 
signal-to-noise  ratio  (dB),  one  interfering  source  at  10°  N  50. 
No  =  200,  Hann  weighting  DOAr  ==  2° 
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Table  1 

Input  parameters  for  the  simulated  tests 


Parameter 

Symbol 

Value 

Fig.  no. 

number  of  statistical  draws 

Nr, 

200 

2-9 

number  of  hydrophones 

L 

16 

2-9 

spacing/ wavelength 

d/\ 

0.5 

2-9 

noise  power  (dB) 

<rl 

0 

2-9 

time-bandwidth  product 

N 

f  50 

l  200 

2-4, 6-8 

5,9 

array  shading  (CVB) 

f  Hann 
l  uniform 

2, 3, 5-9 

4 

Example  1 

number  of  signals 

I 

2 

2-5 

signal- to- noise  ratio  (dB)  j 

,  a\lal 

6  [-20,13.5] 

30 

2-5 

2-5 

signal  location  (deg)  j 

DOAr 

doa2 

6  [-20,9] 

10 

2-5 

2-5 

Example  2 

numv  . :  of  signals 

1 

3 

6-9 

f  *\!*\ 

5,10,15 

6-9 

signal-to-noise  ratio  (dB)  < 

30 

6-9 

l  *11*1 

30 

6-9 

(  DOA, 

€  [DOAj, DOA,] 

6-9 

signal  location  (deg)  ( 

|  DOAj 

1  DOAj 

-10 

-5,0,5 

6-9 

6-8 

l  DOA, 

-5 

9 

Table  2 

Empirical  probability  of  false  alarm,  Example  1 


Algorithm 

Pfa  (%) 

C'VB 

1.1 

AIC 

2.7 

MDL 

0.0 

GLRT 

0.1 
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-2°  (d).  For  the  same  bearing  variation  the  performance  of  the  high-resolution 
algorithms  is  relatively  constant,  unless  for  case  (d)  where  a  slight  decrease  of  about 
2.5  dB  can  be  noted. 


In  Fig.  3  the  improvement  obtained  by  the  high-resolution  algorithms  is  clearly  seen 
in  terms  of  source  location.  As  expected,  the  CVB  performance  is  always  better  than 
that  of  the  other  algorithms  when  the  source  is  far  away  from  the  interfering  source. 
However,  for  a  low-level  source  with  SNRin  =  -6  dB  (a)  the  sidelobe  decreases 
the  CVB  performance  19°  away  from  the  interference  location.  Figure  3b  shows  a 
similar  behaviour  for  an  SNRj„  =  0  dB:  in  this  figure  the  influence  of  the  mainlobe 
appears  for  |DOAt|  <12°. 

In  Fig.  4  the  conventional  beamformer  performance  is  shown  with  uniform  shading 
instead  of  Hann  shading.  Only  one  signal-to-noise  rat:  is  show"-  SNP;n  =  5  dP. 
As  expected,  the  detection  performance  curve  of  the  conventional  beamformer  has 
a  spikey  aspect  which  is  directly  related  with  the  high-level  source  beampattern: 
the  maxima  of  detection  corresponding  to  the  nulls  and  the  minima  of  detection  to 
the  peaks  of  the  beampattern.  The  main  difference,  compared  to  the  Hann  shading 
results,  is  the  fact  that  large  losses  of  detection  can  be  obtained  even  for  directions 
‘far  away’  from  the  interference  direction. 

In  Fig.  5  the  ‘time-bandwidth  product’  was  increased  to  IV  =  200.  The  conventional 
beamformer  performance  (with  Hann  shading)  is  approximately  the  same  while  the 
high-resolutior  .igorithins  behave  significantly  better. 


4.2.  EXAMPLE  2:  TWO  HIGH-LEVEL  INTERFERING  SOURCES 


The  results  of  this  test  are  shown  in  Figs.  6  to  9.  In  these  figures  the  percentage  of 
detection  probability  is  given  vs  the  relative  location  of  the  low-level  source  to  the 
fixed  high-level  source  normalized  by  the  separation  of  the  two  high-level  sources. 
The  fixed  high-level  source  is  located  at  DOAn  =  -  10°.  The  separation  of  the  two 
high-level  sources,  ADOA, ,  is  variable  and  takes  the  following  values.  5°  (Figs.  6 
and  9),  10°  (Fig.  7)  and  15°  (Fig.  8).  Thus,  the  normalized  DOAn  is  given  by 


DOA„ 


DOA,  -  DOAn 

ADOAi 


where  DOA,  is  the  location  of  the  signal  under  detection. 

An  empirical  estimation  of  the  false  alarm  rate  for  the  f nir  algorithms  gave  the 
results  summarized  in  Table  3.  One  can  note  that  the  overall  Ppa  decreased  in  this 
test  and  is  now  situated  around  0.7%. 
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Fig.  3a.  Estimated  probability  of  detection  (%)  vs  the  relative 
location  to  the  interfering  source  at  -10°.  .V  r  50,  No  =  '200, 
Han n  weighting  SNRjn  =  -6  dB 


Fig  36.  Estimated  probability  of  detection  f%!  >’•«  the  relative 
location  to  the  interfering  source  at.  10'  \  ">0.  /Vp  200. 

Han u  weighting  SNRj„  —  0  dB 
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Fig.  4.  Estimated  probability  of  detection  (%)  vs  tile  relative 
location  to  the  interfering  source  at  -10°,  N  -  50,  ND  =  200, 
uniform  weighting  SNRj„  =  5  dB. 


Comments  Figures  6  through  9  show  a  similar  behaviour:  the  high-resolution  al 
gorithnis  could  achieve  some  detections  (and  in  some  cases  a  good  detection  rate) 
whereas  the  conventional  beamformer  did  none.  As  expected,  from  the  results 
of  Example  1,  higher  detection  probabilities  where  obtained  for  larger  SNR’s,  for 
larger  separations  betwwen  the  two  high-level  sources  and  with  maxima  at  the  half¬ 
separation  of  the  interfering  sources.  Also,  as  before,  a  larger  number  of  averages 
improves  the  results  of  the  high-resolution  techniques  (Fig.  9). 


Table  3 

Empirical  probability  of  false  alarm,  Example  2 


Algorithm 

Pfa  (%) 

CVB 

0.6 

A  If 

1.8 

MDL 

0.0 

GLRT 

0.6 
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Fig.  5a.  Estimated  probability  of  detection  (%)  vs  the  relative 
location  to  the  interfering  source  at  -10°,  N  =  200,  Nd  =  200, 
Hann  weighting  SNRln  =  9  dB. 


Fig.  5k  Estimated  probability  of  detection  {%)  vs  the  relative 
location  to  the  interfering  source  at  -10°,  N  =  200,  No  =  200, 
Hann  weighting  SNRjn  —  -6  dB. 
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Fig.  6a.  Estimated  probability  of  detection  {%)  vs  the  normalized 
relative  location  to  the  two  interfering  sources  at  -10°  and  -5°, 
N  =  50,  /Vq  =  200,  Haim  weighting  SNRirt  =  5  dB. 


Fig.  6b  Estimated  probability  of  detection  (%)  t'5  the  normalized 
relative  location  to  the  two  interfering  sources  at  -10°  and  5°, 
N  50,  JV|)  200.  Hann  weighting  SNR,„  10  dB 


21 


SaCLANTCEN  SM-2U8 


Fig  7a.  Estimated  probability  of  detection  (%)  vs  the  normalized 
relative  location  to  the  two  rfering  sources  at  -10°  and  0°, 
N  =  50,  No  ~  200,  Hann  weighting  SNRj„  -  -  5  dB 


Fig.  7b  Estimated  probability  of  detection  (%)  the  normalized 
relative  location  to  the  two  interfering  soiirces  at  -10°  and  0°, 
N  -  50,  No  ~  200,  Hann  weighting  SNRi„  =  0  dB 
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Fig.  8a.  Estimated  probability  of  detection  (%)  vs  the  normalized 
relative  location  to  the  two  interfering  sources  at  -10°  and  5°, 
N  =  50,  ATd  —  200,  Hann  weighting  SNRi„  =  -10  dB 


Fig.  6b,  Estimated  probability  of  detection  (%)  the  normalized 
relative  location  to  the  two  interfering  sources  at  -10°  and  5% 
N  —  50,  Nq  ~  200,  Hann  weighting  SNRjt,  -  -5  dB. 
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Fig.  9a.  Estimated  probability  of  detection  (%)  vs  the  normalised 
relative  location  to  the  two  interfering  sources  at  -10°  and  -5°, 
N  =  200,  No  =  200,  Hann  weighting  SNRi„  =  0  dB. 


Fig.  9b.  Estimated  probability  of  detection  (%)  vs  the  normalized 
relative  location  to  the  two  interfering  sources  at  —10°  and  —5°, 
N  =  200,  No  =  200,  Hann  weighting  SNRi„  =  5  dB 
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5.  Conclusion 

The  problem  of  detecting  a  low-level  source  in  the  presence  of  high-level  interfering 
sources  has  been  studied  in  this  report  using  simulated  data.  The  performance  of 
the  conventional  beamformer  has  been  compared  with  three  high-resolution  beam- 
formers  based  on  eigenstructure  decomposition  techniques.  From  the  simulations 
it  has  been  possible  to  estimate  the  probability  of  detection  of  the  low-level  source 
keeping  the  probability  of  false  alarm  approximately  constant.  As  expected,  the 
simulation  shows  that  the  conventional  beamformer  performance,  although  optimal 
in  a  single-source  scenario,  is  greatly  reduced  when  attempting  to  detect  a  low-level 
source  near  to  high-level  interfering  sources. 

The  detection  performance  of  the  high-resolution  algorithms  is  significantly  better 
than  that  of  the  conventional  beamformer  for  a  low-level  source  located  within  the 
mainlobe  of  a  high-level  interfering  source.  When  more  than  one  high-level  source  is 
present  the  performance  improvement  can  be  extremely  high.  In  the  case  of  a  very 
low-level  source  the  high-resolution  techniques  can  achieve  a  detection  improvement 
even  in  the  sidelobe  region.  The  detection  behaviour  of  the  high-resolution  tech¬ 
niques  is  relatively  constant  in  bearing  except  for  the  interference  direction;  this 
behaviour  is  improved  by  a  longer  observation  time. 

Among  the  tnree  high- resolution  detection  algorithms,  AIC  always  displayed  the 
best  performance  -  together  with  the  highest  false  alarm  rate  -  and  MDL  displayed 
the  best  results  -  together  with  the  lowest  false  alarm  rate.  This  is  in  agreement 
with  the  asymptotical  consistency  remarked  by  [12,17]. 

In  conclusion,  the  results  suggest  that  the  high-resolution  algorithms  can  achieve 
a  real  improvement  in  detecting  low-level  targets  in  ship-induced  noise  fields.  A 
disadvantage  of  the  high-resolution  algorithms  is  their  dependence  on  a  knowledge 
of  the  noise  field  Deviations  from  the  assumed  noise  field,  which  may  occur  in 
real  data,  can  lead  to  a  degradation  in  detection  performance.  In  a  companion 
report  [7]  we  present  real  data-processing  results  which  show  the  robustness  of  the 
high- resolution  detection  algorithms. 
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